Three-particle collisions in quantum wires: 
Corrections to thermopower and conductance 



O 
O 

(N 



c3 



i 

-a 
c 

o 
o 



> 

(N 

o 
m 
o 
l> 
o 

-I— > 



-a 
c 

o 
o 



Anders Mathias Lunde, 1 ' 2 Karsten Flensberg, 1 and Leonid I. Glazman 2 

1 Nano- Science Center, Niels Bohr Institute, University of Copenhagen, DK-2100 Copenhagen, Denmark 

2 William I. Fine Theoretical Physics Institute, University of Minnesota, Minneapolis, MN 55455, USA 

(Dated: February 4, 2008) 

We consider the effect of electron-electron interaction on the electron transport through a finite 
length single-mode quantum wire with reflectionless contacts. The two-particle scattering events 
cannot alter the electric current and therefore we study the effect of three-particle collisions. Within 
the Boltzmann equation framework, we calculate corrections to the thermopower and conductance 
to the leading order in the interaction and in the length of wire L. We check explicitly that the 
three-particle collision rate is identically zero in the case of several integrable interaction potentials. 
In the general (non-integrable) case, we find a positive contribution to the thermopower to leading 
order in L. The processes giving rise to the correction involve electron states deep in the Fermi sea. 
Therefore the correction follows an activation law with the characteristic energy of the order of the 
Fermi energy for the electrons in the wire. 



I. INTRODUCTION 

Short clean one-dimensional (ID) mesoscopic wires, of- 
ten referred to as quantum point contacts, show con- 
ductance quantization^ 2 - as a function of the channel 
width. The quantization is well described by the the- 
ory of adiabatic propagation of free electrons.— For non- 
interacting particles, conductance quantization should 
occur in longer channels too, as long as there is no 
backscattering off inhomogeneities within the channel. 

A lot is known about the role of electron-electron in- 
teraction of ID channels. Electron-electron repulsion in 
a wire enhances dramatically the reflection coefficient, 
making it energy-dependent^ However, interaction be- 
tween electrons does not alter the quantization (in units 
of 2e 2 /h) of an ideal channel conductance in the limit 
of zero temperature.^ What is still an open question 
is whether there are other manifestations of interactions 
due to inelastic processes, which influence the transport 
properties. 

In the absence of interactions, left- and right-moving 
particles in a wire are at equilibrium with the reservoirs 
they originate from. If a bias is applied between the 
reservoirs, then these equilibria differ from each other, 
giving rise to a particular form of the non-equilibrium 
distribution inside the channel. On the other hand, in 
a long ideal channel and in the presence of interactions 
one may expect equilibration to occur between the left- 
and right-movers into a single distribution characterized 
by an equilibrium with respect to a reference frame mov- 
ing with some drift velocity. Interestingly, in a model 
with momentum-independent electron velocity for left- 
and right-movers (as it is the case in Tomonaga-Luttinger 
model) there is no difference between the two distribu- 
tions. Effects originating from the particle-hole asymme- 
try, however, may discriminate between the two. Ther- 
mopower and Coulomb dragZ^ are examples of such 
effects. 

At present, little is known about equilibration in a 
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FIG. 1: A schematic picture of two metallic gates deplet- 
ing the underlaying two dimensional electron gas and thereby 
forming a short ID quantum wire of length L. This fabri- 
cation method has the advantage of producing reflectionless 
contacts to the leads,— so that the boundary conditions of the 
distribution function is given by the Fermi function of the 
reservoirs. We define the thermopower as S = V/AT| , 
i.e. the voltage V required to counteract a current due to the 
temperature difference AT. 



ID electron system. In higher dimensions the electron- 
electron interaction provides the most effective relaxation 
mechanism at low temperatures and therefore we include 
this relaxation mechanism as the first approach. How- 
ever, in ID pair collisions cannot change the distribu- 
tion function for quadratic dispersion, since the momen- 
tum and energy conservation^ laws result in either zero 
momentum exchange or an interchange of the two mo- 
menta^ In either case the distribution function remains 
the same. Thus, the leading equilibration mechanism is 
due to three-particle collisions, which we study in this 
paper fi 2 - 

We investigate here the effects of three-particle col- 
lisions in reasonably short wires (see Fig. Q]), where 
electron-electron scattering can be considered perturba- 
tively. As measurable quantities, we evaluate the temper- 
ature dependence of the thermopower and conductance. 
Note that for more than one mode pair collisions become 
important for certain fillings 

The paper is organized as follows: First we review the 
non-interacting limit of thermopower and give a qualita- 
tive explanation of the effects due to three-particle col- 
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lisions. Then we describe how to include the electron 
interactions using the Boltzmann equation. Next we cal- 
culate the main ingredient for our perturbation theory, 
namely the three-particle matrix element and scattering 
rate using a T-matrix expansion. We note several inter- 
esting properties of this scattering rate. Finally, we de- 
rive the conductance and thermopower corrections and 
discuss the deviation from the so-called Mott formula. 
Furthermore, some technical details are put in two ap- 
pendices, and in Appendix [A] we show that the number 
of left and right movers have to change in a scattering 
event for the current to change. 



A. Thermopower in the non-interacting limit 




For a wire without interactions the distribution func- 
tion is determined solely by the electron reservoirs 



(0) 



/°( £fc - ML ,T L )^/ L °( £fc ) forfc>0, 
/°(e fc -M R ,T R ) = /°( £fc ) forfc<0, 



(1) 



where e k is the dispersion relation for momentum k and 
spin a (suppressed in the notation), and f°(e — fx,T) = 
[l + cxp(( £ — ^)/fc B T)]~ 1 is the Fermi function with /Ul/r 
and T L / R denoting the chemical potential and tempera- 
ture of the left /right contact, respectively (see Fig. [J). 
The electric current for low temperature T <C Tp and in 
linear response to the applied bias V and temperature 
difference AT -C T then follows as (e > 0)^ 



I (o )= (_Aj2v k (f°(e k )-f°(e k )) 



(2) 



crfc>0 



- ~V (l - 6 - Tf/T ) + f k B AT^e-^ T . (3) 

From this the well-known leading-order results for con- 
ductance 



G (0) = ^- (l - e- T *' T ) 
and for thermopower 



g(o) _ c -t f /t 
e T 



(4) 



(5) 



for a fully open channel are obtained. Here Tp = £ F /k B 
is the Fermi temperature. 



FIG. 2: (a) The dominant three-particle scattering process 
at low temperature in a single energy band, (b) The three- 
particle scattering process perturbing the initial distributions 
shown with warm left moving electrons (k < 0) and cold right 
moving electrons (k > 0). Due to the temperature difference 
of the initial distributions the scattering process creating left 
movers dominate compared to the opposite scattering and 
therefore it gives a positive correction to the thermopower. 



current the number of left and right moving electrons 
need to change, since it is the number of electrons going 
through a mesoscopic structure that determines the cur- 
rent and not their velocity (see the Appendix [A| . Sec- 
ondly, we find the dominant scattering process at low 
temperature to only involve a single electron changing 
direction. This occurs near the bottom of the band as 
pictured on Fig. [2](a) . For the initial electronic distribu- 
tion the left moving electrons have a higher temperature 
then the right moving ones, which favors scattering into 
the warmer distribution as seen on Fig. EJb). This thus 
creates more left moving electrons and thereby increases 
the particle current towards the colder reservoir, i.e. in- 
creasing the thermopower. 

Another important point is that the thermopower and 
conductance corrections are exponential in temperature, 
i.e. proportional to exp(— Tp/T). This is a direct conse- 
quence of the dominant three-particle scattering process 
requiring an empty state near the bottom of the band. 
We find the form of the thermopower correction at low 
temperatures due to the three particle scattering to be 
given by 



B. Main results and a simple picture of the effect 
of the three- particle scattering 

One of the main results of this paper is that the three- 
particle collisions give a positive contribution to the ther- 
mopower, i.e. the current due to a temperature differ- 
ence is increased by the three-particle scattering. This 
can be explained in simple terms. Firstly, to change the 
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T 



cxp(-T F /T) > 0, 



(6) 



where V is the electron-electron interaction strength and 
Tp the Fermi temperature. This is found perturbatively 
in the short wire limit. The long- wire limit remains an 
open question, and we expect that the length dependence 
of thermopower saturates once L exceeds some relaxation 
length (which increases for decreasing temperature). 
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In contrast, the conductance correction is negative. To 
understand this, note that the chemical potential of the 
initial distribution is higher for the right moving electrons 
then the left moving ones. This favors scattering into 
the left moving branch (still with the process shown in 
Fig.[2ta)) for non-zero temperature and thereby decreas- 
ing the current. The form of the conductance correction 
is similar to the thermopower correction: 



where fk{$) is the distribution function at a space point 
x between zero and L (see Fig. [TJ, v k = jid k e k is the 
velocity and T kx [/] is the three-body electronic collision 
integral, i.e. no impurity or interface roughness effects are 
included here. We include the voltage and temperature 
difference in the boundary conditions of the reflectionless 
contacts,— i.e. 



G" 



-L\V\ A 



T 



cxp(-T F /T) < 0. 



(7) 



II. CURRENT CALCULATION IN THE 
BOLTZMANN EQUATION FORMALISM 

A. Effect of interactions on the current 

To model the current through a short ID quantum wire 
including perturbatively the three particle interactions, 
we use the Boltzmann equation 



Vkd x fk{x) = Xkx[f], 



(8) 



fk(x 
fk(x-- 



= 0) 
L) 



for 
for 



k > 0, 
k < 0, 



(9a) 
(9b) 



and therefore omit the term kd k fk{x) in the Boltzmann 
equation, allowed in the linear response regime.— A sim- 
ilar method has been used to investigate electron-phonon 
interactions in short quantum wires,— quantum Hall ef- 
fect in quantum wiresii and ballistic Coulomb dragi^ 

The three-particle collision integral is assumed to be 
local in space and is given by 



i k Af] = - 



w 123 ;i>2> 3 > [/i/ 2 / 3 (i - - /a0(i - h>) - h'hh{i - /0(i - h){i - h)] , (io) 



a ^ a 3 k 2 k 3 
<7 lt cr 2 ,a 3l ki , k2 , k3 



where the quantum numbers are primed/unprimed af- 
ter/before the scattering event, /; = fki(%) and the 
scattering rate Wi23 ; i'2'3' is found in the next section. 
Without interactions (Wi23;i'2'3' = 0) the solution of the 
Boltzmann equation is simply given by /(°) in Eq. (fT|). 
When interactions are included it becomes a very difficult 
task to solve the Boltzmann equation to all orders in the 
interaction. However, for a short wire the interactions 
only have short time to change the distribution function 
away from the initial distribution and therefore we 
expand the distribution function in orders of W r i23 ; i'2'3 / 
as 



f = f (o) + fW + 



(i) 



(11) 



To find to the first order in W, we insert the ex- 
pansion of / in the Boltzmann equation and realize that 
only f(°> is necessary in the collision integral. Since 
1kx[f^\ = 2fc[/(°)] is independent of x, we find that 



for k > 0, 
I k [f m ] for fc<0, 



(12a) 
(12b) 



using the boundary conditions Eq. (J9j) . Therefore the 



current to the first order in W is 

/ = j(°) + e Y, i k [f (0) ] = i (0) + r a \ 



(13) 



erfc<0 



where 1^ is the non-interacting (Landauer) part of the 
current from Eq. (f2]) and / lnt is the part due to interac- 
tions. 



B. The linear response limit 

The form of the interacting part of the current is now 
known and the next step is therefore to evaluate it to 
linear response to V and AT to obtain the thermopower 
and conductance corrections. To this end, we define ■0^.°^ 
via 



/( 0) S / ( £fc ) + / ( £fc )(l-/°(e fc ))Vr 



CO) 



(14) 



where f°{sk) is the Fermi function with temperature T 

and Fermi level e F . It turns out that ip^' is proportional 
to either V or AT. This is seen by using the identity 



-k B Td E f(e k ) = f(e k )(l-f°(e k )) : 



(15) 
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so we can identify tp^ by expanding the non-interacting have used the relation 
distribution function (see Eq. ([TJ and Fig. Q]) 

fl(e k ) ~ /°(e fe ) + [-a e /°(e fc )] eV, (16a) 



AT 

/°( £fc ) ~ /°( e *) + [-a e /°(e fc )] (e - £ F ) — , (16b) 



( eV 

~ \ efc-£p AT 
I k„T T ' 



for k > 
for fc < 



(17) 



Therefore to get T nt in linear response to V and AT, 
we linearized the collision integral X fc [/( )] Eq. (fit))) with 

respect to and insert it into I mt Eq. (|13|1 to obtain 



' int = HoE E 



^123:l'2'3' 



where we defined 



^(o)_^(o)_ V) (o) 



(18) 



Ai23 ; l'2'3' — 

WW 2 '3' A /2°/ 3 (1 - fv) (1 - ./2 ' )(1 - /&), (19) 



using the shorthand notation ipf^ = and /P = 

/ ( £ fc.)- To linearize the collision integral and thereby 
the correction to the current due to interactions T nt , we 



do) 



/l /2 /3 (l-/l ')(l-/2')(l-/ 3 ") = 

/l '/2 '/3 '(l-/l )(l-/2°)(l-/ 3 ), (20) 

valid at e\ + e 2 + £3 = £i' + £2' + £3' • 

Since ipi°^ ls different for positive and negative fcj, we 
need to divide the summation in T nt Eq. I|18p into posi- 
tive and negative k sums, which gives 2 5 = 32 terms. For 
this purpose, we introduce the notation 

E (■)=£(■), £(•) = £(•), (21) 



fci<0,fc 2 >0,fc 3 <0 
fc 1 />0,fc 2 />0,fe 3 /<0 



++- 



spin 



and similarly for other combinations of the summation 
intervals. The 32 terms can be simplified to only three 
terms using energy conservation and symmetry proper- 
ties of Ai23 ; i'2'3' in Eq. (fl9|) under interchange of indices. 
There are pairwise exchanges of indices Ai23 ; i'2'3' = 
A2i 3; i'2'3' = Ai23 ; i'3'2' etc. and interchanges between 
primed and unprimed indices, A 12 3 ; i'2' 3 ' = A 1 / 2 '3' : i23, 
using Eq. (|20|l and the fact that Wi23;i'2'3' contains a 
matrix element squared. This leads to six terms. Fur- 
thermore, Ai23 ; i'2'3' is invariant under fcj — > — ki for all 
i = 1, 2, 3, 1', 2', 3' simultaneously due to time reversal 
symmetry, also seen explicitly from the form of Wi23;i'2'3' 
(derived below) . An example of how the simplifications 
occurs can be seen in Eq. l|A6| . Thus we obtain the more 
compact result 



T nt = 2(-e)E E A i23;i'2', 



spin — I — h 
+++ 



AT eV 

(e 1 - £ F ) 



k B T2 



k B T 



4 (- e )£ E Ai23 ^' 



spin h 

++ + 



3 (- e )E E Ai ^!' 

spin -H — 



2'3' 



2'3' 



AT . . 2eV 

{(£l -£p) + (£2 -£ F )} - 



k B T* 



k B T 



AT eV 

{-(£3 - £ F ) + (£2' - £ F ) + (£3' - £ F )} 



k B Ti 



(22) 



where the definition of ip\ in Eq. fT7|) was inserted. An 
important point is that the number of positive/negative 
wave- vector intervals is not the same before and after the 
scattering. Therefore, we note that only scattering events 
that change the number of left and right moving electrons 
contributes to the interaction correction to the current. 
The origin of this is the cancellation of the velocity in the 
definition of the current and in the distribution functions 

This cancellation thus leads to an expression for the 
interaction correction to the current in Eq. ([22| where 



all the in-going and out-going momenta enter on equal 
footing. In the Appendix [Aj we show that this is valid 
to all orders in perturbation theory. Due to this prop- 
erty and momentum conservation, there are no processes 
that alter the current possible near the Fermi level. Con- 
sequently, states far away from the Fermi level have to 
be involved in the scattering, which, as we will see, leads 
to a suppression of T nt by a factor exp(— T F /T). The 
distribution function on the other hand, can be changed 
by scattering processes near the Fermi level. 

To identify the important processes we find in the next 
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section the scattering rate Wi23 ; i'2'3'- 



III. THREE-PARTICLE SCATTERING RATE 

The three particle scattering rate W / i23;i'2'3' is calcu- 
lated using the generalized Fermi Golden rule inserting 
the T-matrix, T = V + VGoT, iterated to second order 
in the interaction V to get the three particle interaction 
amplitude, i.e., 



Wl23; 



2- 



1'2'3' 



|<l , 2'3'|VG ^|123) c | 2 <y( J B i - Ef), (23) 



where E i = e 1 + £ 2 + e 3 m the initial energy, E* = 
Ey + e-2' + £3' the final energy, Go is the resolvent op- 
erator (or free Green function), j is short hand for kj, 
and the subscript V means connected in the sense that 
the scattering process cannot be effectively a two parti- 
cle process, where one of the incoming particles does not 
participate in the scattering. Explicitly Go and V are 
given by 



Go 



1 



Ei-Ho + irj' 



(r) 0+) 



(24) 



V ~ 2L E ^ V « C k 1 +q<? 1 C i 2 -qa 2 C k 2 * 2 c k 1 a 1 - ( 2 5) 
fe 1 fe 2 ij o x o 2 

Here Hq is the unperturbed Hamiltonian (i.e. kinetic en- 
ergy with some dispersion), V q the Fourier-transformed 
interaction potential, and c k(J (4a-) is the annihilation 
(creation) operator. To calculate the matrix element 



(l'2'3'|VGoV|123) c , we write the initial and final states 
as 



|i23) = 4 i(Ti 4 2(72 4 3(r3 |o), 



(26) 

|l / 2 , 3 / ) = 4 1 ,. 1 ,4 2 ,. 2 ,4 3 , CT JO), (27) 

where |0) is the empty state. Using the anti-commutator 
algebra {c^c^} = S^j, we obtain 



G y|123)=^^K, 



(28) 



E 



(abc)eP(123) 

where we introduced 



sgn (abc) 

£ba(q) 



4 a +ga Q 4 -go- 4 c cr c |0)> 



£ba{q) = £b + £a - £b-q ~ £ a +q + if] 

= — q(h - k a - q) + irj 



(29) 



(the last equality is only valid for a quadratic dispersion) , 
and where the set of permutations is given by P(123) = 
{(123)+, (231)+, (312)+, (132)-, (321)-, (213)-}. Here 
the signs of the permutation, sgn(a6c), are shown as su- 
perscripts. 

In order to exclude the effectively two-particle pro- 
cesses when multiplying Eq. I|28p by (1'2'3'|V from the 
left, k c (c — 1,2,3) needs to be different from kji 
(j = 1,2,3). The result is 



(1'2'3'|VG V|123) C 



E 



(2L)2 ^ 

v \abc)eP(123) (a'6'c')eP(l'2'3') 



Vi i \ / lit i\ Va'—aVc' — c 3a+b+c,a'+b'+c' err 
sgn(a6c)sgn(a be )- — — : rrz"a n , ,<?„o<T h , ,<j k o a 



£b 



Eb+c 



"I 



(30) 



where V q = V q + V- q is the symmetrized interaction. The matrix element consist of 36 terms and the scattering rate 
thus has 36 2 = 1296 terms. To obtain this result we did not use energy conservation. For a quadratic dispersion, the 
denominator is only zero if we have an effective pair collision or if the momentum transfer is zero, as seen from the 
expression £ba{q) = ^q(kb — k a — q) + irj. A picture of the matrix element is found in Fig. EJb) , where the exchange 
processes (including the sign) are visualized as different ways to connect two interaction lines and an intermediate 
propagation (Go) seen on Fig. [3ja). The inclusion of the Fermi statistics makes a substantial difference for the 
properties of the scattering rate as compared to the case described in Ref. which is obtained by setting all 

sgn(- ••)=+!• 

We can rewrite the matrix element Eq. (pO)) in a more transparent way in terms of quantum mechanical exchange 
symmetry. First, we introduce the following combination of three-particle scattering amplitudes 



V(] 



4L 2 

Vl'-lV2'-2 



Vl'-lVy-\ 



V?,' -3V2' -2 



£3 + £2 — £3' — £2+3-3' £l + £3 ~ £l' — £3+1-1' £2 + £l — £2' — £l+2-2' 



T -2 + £3 — £2' — £3+2-2' 



+ 



V21-2V31-3 



r i + £2 — £v — £2+1-1' £3 + £1 — £3' — £1+3-3' 
I 



(31) 



and after some rewriting, we then obtain 
(l'2'3'|l/G y|123) c = 6 kl+k2+k3tkll+k2l+ksl (32) 
V(ll', 22', 33') + V(12', 23', 31') + V(13', 21', 32') 

V(ll', 23', 32') - V(13', 22', 31') - V(12', 21', 33') 



We interpret this result in a way similar to a two-particle 




FIG. 3: A visualization of the connected three-particle scattering matrix element Eq. I|30p. where three particles interchange 
their momenta and energy. This matrix element enters the scattering rate via the generalized Fermi Golden rule Eq. lf23|) . (a) 
The basic three-particle interaction consisting of two interaction lines and a free propagation, see Eq. lf30|l . (b) Picture of the 
exchange processes times the basic interaction needed to form the matrix element Eq. (|30[) . 



matrix element, 



<1'2'|V|12) = 



^4 1 ,-fc 1 ( ^ 1 ,ov < ^ 2 ,er 2 



Vu 



, (33) 



which contains a direct (first term) and a exchange term 
(where V <-► 2'). 

In the three-particle case V(ll', 22', 33') is the direct 
term and one can make five exchange processes (instead 
of one) by exchanging the three final states 1', 2' and 
3'. This gives Eq. (021). The sign in front of each ¥(•■•) 
is determined by the number of exchanges made, e.g. in 
V(ll', 23', 32') a single exchange, 2' <-> 3', gives a minus 
(-1) 1 whereas for V(12', 23', 31') two exchanges (1' <-> 3' 
followed by 3' «-> 2') gives a positive sign (— l) 2 . Fur- 
thermore, the arguments in V(ll', 22', 33') are ordered 
in three pairs such that the differences between the el- 
ements in each pair are the only arguments of the in- 
teraction potential, see Eq. IpTjl . This is useful when 
constructing approximations having a specific scattering 
process in mind. 

How the matrix element was rewritten into the form of 
Eq. ([321 can also be described in terms of the drawings 
of Fig. [3 The direct term V(ll', 22', 33') is the sum of 
the six terms having mirror-symmetric exchanges before 
and after the scattering. The other terms in Eq . (|32)1 then 
can be obtained by suitable changes of outgoing lines. 



Remarkably, for some two-body potentials, scattering of 
the particles of an JV-body system is exactly equivalent 
to a sequence of two-body collisions. Such "special" po- 
tentials were studied in the context of integrable quan- 
tum many-body problems <^ We recall now that for a 
quadratic band a pair collision does not change the mo- 
menta of the incoming particles or simply permutes the 
two momenta. Therefore, three-particle scattering for a 
the integrable potentials may result only in permutations 
within the group of three momenta of the colliding par- 
ticles; all other three-particle scattering amplitudes must 
be zero for such potentials. In the context of this work it 
means that even three-particle (or higher-order) collisions 
would not bring electron equilibration for such types of 
electron-electron interaction. 



In this section we check that the three-particle scatter- 
ing amplitudes are indeed zero for two special potentials. 



Zero three-particle scattering rate for 
integrable models 



Point-like interaction 



The expressions we obtain for the three-particle scat- 
tering rates Eq. lj23|) are quite cumbersome. Nevertheless, 
the obtained results allow for some consistency checks. 



In the case of contact interaction, V„ 



constant 



Vq, and for any kind of electron dispersion relation (i.e., 
not necessarily quadratic), we find by using the energy 
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conservation law that 



^|(l'2'3'|^GoF|123) c | 2 

spin 



{2LY 

|Ai21' — Ai22> — A\SV - 
+ \Ai2V 

+ \A 122 
+ \A 12 i 



0~k 1 +k 2 +k 3 ,k 1 , +k 2 , +k 3 , 



1132' 



\A 13V 
\A 12 v 

\Ai3V 

\A 122 < 
\A 132 > 



Ai23> — A131' 
Ai23' — A132' 
A\22' — A 2 3\i 
A\Z2' — ^231' 
^123' — ^231' 
^4l33' — ^231' 
^123' — ^232' 

^133' — ^232' 



133' 



133' 



232' 



232' 



233' 



23:',' 



233' I 



A 



233' 



(34) 



where A abc = (e a + s b - e c - £ a +b-c + 1 . This is 
a major simplification from 36 2 = 1296 to 9 x 4 2 = 144 
terms by performing the spin summation. If furthermore, 
the dispersion is quadratic, Sk oc fc 2 , then we find the (at 
first sight) surprising cancellation 



5^|(l / 2 / 3 / |VG V|123) I 



x S(ei + e 2 + e 3 - e v - e 2 > - £3') = 0. (35) 



This can be seen directly from Eq. (|34j) or by noting that 

V(lo', 26', 3c')5(ei + e 2 + e 3 - e a > - e v - e d ) = 0, (36) 

for a quadratic dispersion and constant interaction for 
(a'b'c') e P(l'2'3') ) i.e. each term of Eq. (02]) is zero. 
V(ll', 22', 33') cancels in such a way, that the three first 
terms of Eq. (|3"Tj) cancel each other (the even permuta- 
tions of (123) combined with the same primed permuta- 
tion) and the three last terms cancel each other (the odd 
permutations of (123) combined with the same primed 
permutation) . 

In fact, the described above cancellation is in agree- 
ment with the general factorization results for the S- 
matrix of one-dimensional TV-body problem with 6- 
function interaction in real space.— In this context, it is 
crucial that the particles have quadratic dispersion rela- 
tion; if we use e.g. £k oc k 4 , then the cancellation does not 
occur. Notice also that the cancellation we demonstrate 
is not a trivial zero. Indeed, the underlying two-particle 
amplitudes Eq. (|33|) are finite for a ^-independent po- 
tential, if one includes spins. (For spinless fermions and 
contact interaction the matrix element would be zero be- 
cause the direct and the exchange terms cancel in accor- 
dance with the Pauli principle.) 



2. The V q = Vo(l — q /go) interaction 

We checked also that the energy conserving part of the 
matrix element (l'2'3'|VG ! o^|123) c in the case of spin- 



less fermions, quadratic dispersion, and the Fourier trans- 
formed interaction potential of the form 



V„ = V [1 



4 



(37) 



becomes equal to zero. This is also possible to ex- 
pect because of the relation of the potential Eq. lf37|) 
to the integrable ID bosonic Lieb-Liniger model 2 ^. In- 
deed, the bosonic model with contact interaction poten- 
tial oc gB${xi — x 2 ) may be exactly mapped^ 2 - onto the 
spinless fermionic model with interaction Vf{x\ — x 2 ) oc 
— (1/9b)8"(xi — x 2 ). The integrability of the bosonic 
model guarantees the integrability of the corresponding 
Fermionic one. Adding a contact interaction to Vf does 
no harm, as we are considering spinless Fermions. Fi- 
nally, Fourier transformation takes us to Eq. J37|). 

We checked that including the spin degree of freedom, 
spoils the remarkable cancellation for a three-particle am- 
plitude. 

In the following Sections, we assume a general case 
interaction potential for which the three-particle scatter- 
ing amplitudes lead to a non-trivial re-distribution of the 
momenta between the particles. 



IV. THERMOPOWER AND CONDUCTANCE 
CORRECTIONS DUE TO THREE PARTICLE 
INTERACTION 

In this section, we go though the main ideas and ap- 
proximations in evaluating the current correction due to 
interactions P nt Eq. ([22)1 to lowest order in the temper- 
ature, T < Tp. We give a more detailed calculation in 
the Appendix IbI 

As noted previously, all three terms in / lnt Eq. ([22)1 are 
exponentially suppressed, since momentum conservation 



ky + k 2 > 



forbids scattering processes near the Fermi level for the 
given combinations of positive and negative k intervals. 
To be more specific, it is the phase space restrictions of 
the Fermi functions that give the exponential suppres- 



sion, i.e. 



/l /2 /3 (l-/l ')(l-/2 ')(l-/3°')«e 



-Tjf/T 



(38) 



We begin by identifying the most important three parti- 
cle scattering process. The three terms in I mt Eq. (|22|) 
are: (i) two right movers backscattering a left mover 
while remaining right movers; (ii) one right mover keep- 
ing its direction, while backscattering two left movers, 
and (iii) a left and a right mover keeping their directions, 
while backscattering the third particle. From now on, we 
will concentrate on the case of Coulomb interaction V q , 
which is the largest for small q, therefore we want to iden- 
tify processes where the initial and final state are close in 
momentum spaced. Further, the process(es) should not 
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require more than one electron in states suppressed ex- 
ponentially by the Fermi functions. One can see that due 
to the constraints stemming from momentum and energy 
conservation, in fact only process (hi) allows both initial 
and final states to be close to each other in momentum 
space and at the same time having only a single expo- 
nentially suppressed factor. The corresponding scatter- 
ing process is of the type shown in Fig.[2|a). Therefore to 
the first order in exp(— Tp/T), we include only the third 
one in Eq. ((22]). This leads to 



123:l / 2 / 3 / 



spin +H 



AT 



{-£3 + £2' + £3' - £ F } - 



eV 



(39) 



Here Ai23 ; i<2'3' expresses the available phase space in 
form of the Fermi functions and the three-particle scat- 
tering rate, see Eq. lfl9|) . 

One essential approximation is that for the scattering 
process depicted in Fig. [2](a) , we may replace the full 
Fermi distribution functions by the exponential tales or 
the low-temperature limit expressions, i.e. 



ft 
fS 



0(k F 
0(k F 



fci)0(fci) 
k 2 )9(k 2 ) 



f°y 



= -(£3-£ F )/fc B T 



1 

l-/ 3 °< 



f°2< 



6{k v 

p( e 2 / 



-*f), 

-e F )/fc B T 



D (e 3 ,-e F )/k B T 



(40a) 
(40b) 
(40c) 



Note that ki, ky and k 2 are all positive. We see that the 
product of the Fermi functions is indeed exponentially 
suppressed, i.e. ex exp(— Tp/T). 

The second essential approximation is that for the scat- 
tering process seen in Fig.EJa) the initial and final states 
differ by a small momentum. Therefore the matrix el- 
ement in the transition rate Wi 23 -y 2 ' 3 ' is dominated by 
the direct term V(ll', 22', 33') in Eq. ([32]), since the five 
exchange terms are suppressed by the Coulomb interac- 
tion |V[ g |~fcJ < I^KfcJ, i-e. 

(l'2'3'|FG y|123) c ~ 

6 kl+k2+k3tkl , +k2 , +k3t V(W, 22', 33'). (41) 



The direct term Eq. (j41[) would be zero for V q — const. 
In the case of quadratic dispersion relation and general 
V q , it vanishes in the limit (kj> — ki) — » and (kf — 
kj) — > for i,j 6 {1,2,3} due to the Pauli principle. 
For a quadratic dispersion and for a general symmetrized 
interaction V q = V q + V- q , the direct term V(ll', 22', 33') 
simplifies to the following expression 



V(ll',22',33') = 



(91 + 93) 



4L 2 H 2 /m 

(91 + 93)^^3 + V qi+q3 (q 3 V qi + 9i V q3 ] 



(h -k 3 + qi)qiq 3 (ki - k 3 - q 3 ) 



(42) 



Next, we give a qualitative explanation for the power 
law in T for the interacting current correction Eq. lj39]) us- 
ing the quadratic dispersion. First, we consider the phase 
space constraint. To do the sum over all k in Eq. f39|) 
we use the momentum and energy conservation and in- 
troduce new variables q\ = ky — k\ and q 3 = k 3 > — fc3, 
i.e., change the summation variables, 



k 1 ,k 2 ,k 3 ,ky,k 2 ',ky -> k 1 ,k 3 ,q 1 ,q 3 . 



(43) 



The energy conservation for a quadratic dispersion gives 
a factor of 1 / |<7i + 53 1 (see e.g. Eq. (|B6P ). For the process 
at hand, the k\ and k 3 are close to the Fermi level and 
each of their sums contribute with a factor of qi and q 3 , 
respectively. The Fermi functions give the exponential 
suppression and a contribution to the phase space in form 
of an exponential tail, i.e. 



/i°/ 2 /3 (l-/°)(l-/2°')(l-/3°)« 



e -T F /T e (e 2 



-£3+£ 3 ')/ fe B T 



(44) 



see Eq. (|40aM40c|) . To get the low temperature result 
for / lnt , Eq. ([391 , we use the method of steepest decent 
to calculate the integral. To this end, we note that the 
exponent e 2 > — e 3 + e 3 > is a function of q\ and q 3 and in 
the limit T/T F — > the most important part is around 
the origin qi = q 3 = 0. Here e 2 i — £3 + £3' vanishes as 
— \hv ¥ {qi + q 2 ) (see Appendix [B] for details). Therefore 
collecting the phase space factors the current correction 
due to three particle interactions Eq. ([39| becomes 



jint K £ -T F /T 



dqi / dq 



9193 



I9i + 93 
x |V(11',22',33')| 2 
"AT Xp I _ qi + q 3 
T T 



-1 - 



L (gi+93)/fc F 



eV 



k B T 



in the limit T <C Tp. Furthermore, it turns out that the 
constraints k 2 > and k 2 i < in the sum Eq. ([39[l only 
leaves phase space close to qi = q 3 for T/Tp — > 0, so we 
can set q 3 = <ji in the integrand and do the integral over 
q 3 , which is oc q\ due to the phase space limits. To lowest 
order in temperature this yields 



7 int oc e~ T ^ T 



dq g 3 e-^ 2 «/ fe F|V(ll',22',33')p 



AT Tp 
~T~~T 



eV 



k B T 



(46) 



From this we conclude that phase space alone (i.e. as- 
suming |V(11',22',33')| 2 to be a constant) gives a tem- 
perature dependence of the form 



CX £ 



T 



AT T F 
~T~ T T 



eV 
k^f 



(phase space only). 

(47) 



where we used energy conservation and introduced q\ 

ky — k\ and q 3 = k 3 i — k 3 . 



However, as we have seen the three-particle interaction 
rate has delicate momentum dependence that needs to 
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be taken into account. Therefore, to calculate the di- 
rect interaction term V(ll', 22', 33') we expand the sym- 
metrized potential V q for small q, as 



V„ = V a 



(48) 



where the parameter go *C k F describes the screening 
due to the metallic gates near the quantum wire and Vq 
is (twice) the q = Fourier transform of the Coulomb 
potential cut-off by the screening. Setting q 3 = qi = q 
into the three-particle scattering rate Eq. (|42|) . we obtain 



(49) 



V(ll',22',33')ocl/ 2 ( ^ 

. qo 



to lowest order in q. Inserting this into the Eq. f46j) 
the final result for the current correction, including both 
phase space factors and the momentum dependent scat- 
tering rate, becomes 



T 8 Vn 



■i r 



AT T F , eV 

— y 



k B T 



(50) 



(Here we noticed that the non-constant three particle 
scattering rate gave rise to four extra powers in tem- 
perature.) The detailed calculation given in Appendix iBl 
yields a prefactor, and the end result is 



8505 _If_ e (Vok F ) 



2048tt 4 
T 



T 

H 
AT 



■(^f) " 



eV 

F 



o 



Qo 
T 
Tf" 



(51) 



it gives rise to four more powers, which results in the T 7 
dependence. 

In the limit of a point-like interaction, qo — ► oo, the 
corrections are zero in agreement with the result of sec- 
tion nnsi 

It is known from the Luttinger liquid theory that in 
the limit of linear spectrum, which corresponds to Tp — > 
oo, the conductance remains finite even if the wire is 
infinitely long (L — ► oo). Therefore it is tempting to 
speculate that the two terms in the square brackets of 
Eq. (|53j) are the first terms of an expansion in A = L/£ eee 
of some function Fc(X) which saturates at a constant 
value in the limit A — > oo. One may also have a similar 
speculation generalizing Eq. ([52]) for the thermopower, 
[...] F s (L/£ eee ). 

As a final remark, we note that the so-called Mott for- 
mula^ relating the thermopower to the low-temperature 
conductance, 



S = -k u T— 



1 dG 



Gde, 



(55) 



is clearly violated by Eqs. f53|) and l|52j) . This violation 
could be expected, because the conventional derivation of 
the Mott formula (for the non-interacting case) assumes 
that the main contribution to the conductance and ther- 
mopower comes from the states around the Fermi level, 
in an energy interval of the order of temperature.— How- 
ever, in the considered case the main contribution to S 
comes from the "deep" states, even in the zeroth order 
with respect to the interaction potential. Correspond- 
ingly there is no surprise that Eqs. |52|) and l|53|) be- 
ing substituted, respectively, in the left- and right-hand 
sides of Eq. ([55|) produce a parametrically large mismatch 
~ Tp/T . 



Combining this result with the zero-order in the inter- 
action terms, see Eqs. (@]) and ((5j), we find for the ther- 
mopower and conductance in the low temperature limit, 



1 



_ k B T F c ~t f /t 
e T 



C -_ e r -T F /T 

h h 



L 



L 



L 



Here we introduced the effective length 
relation^ 



! 8505 (V k F ) 4 
eee 2048^ 3 si 



F) 



(52) 
(53) 
by the 

(54) 



which may be viewed as a mean free path with respect 
to backscattering for a hole near the bottom of the band. 

To recapitulate, the temperature dependence T oc T 7 
in Eq. (|54|) can be understood in the following way: the 
three-particle scattering of a single particle leaves five 
free momenta, and since two are taken by energy and 
momentum conservation this gives T 3 . In addition, the 
interaction, V q , is proportional to q 2 , and when squared 



V. SUMMARY AND DISCUSSION 

We have calculated the leading interaction correction 
to the transport properties of a clean mesoscopic wire 
adiabatically connected to the leads, using perturbation 
theory in the length of the wire. 

For a single-mode wire, the leading interaction correc- 
tions turns out to be given by three-particle scattering 
processes. This is because two-particle processes can- 
not change the current due to momentum and energy 
conservation. To calculate the effect of the three-body 
processes, we have utilized the Boltzmann equation for- 
malism, with three-particle scattering events defining the 
collision integral. We have identified the leading-order 
scattering processes and found that they involve at least 
one state near the bottom of the band, i.e. far from 
the Fermi level. The involvement of such "deep" states 
results in an exponentially small, oc e~ Tp l T , interaction- 
induced correction to thermopower and conductance at 
low temperatures. 

The account for interaction in this paper is performed 
for relatively short wires, where perturbation theory in 
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the interaction or equivalently in the wire length is valid. 
For longer wires one needs to find the distribution func- 
tion by treating the collision integral in the Boltzmann 
equation non-perturbatively. It is not clear whether the 
relaxation of the distribution function would instead yield 
non-exponential corrections to the transport coefficients 
for longer wires. However, since the scattering processes 
that contribute to the current must involve a particle that 
changes direction (which is proven in Appendix [A} , one 
might speculate that the exponential suppression is valid 
for all lengths, as long as electron-electron scattering is 
the only active relaxation mechanism. 

The question of what the relaxed distribution func- 
tion looks like for a mesoscopic wire is an interesting and 
unsolved problem. Here we have only given a partial an- 
swer for the leading contributions for a short wire, i.e. 
to lowest order in the interaction. Further studies should 
involve a self-consistent determination of the distribution 
function. 

Since thermopower is sensitive to the electron distri- 
bution function, it might be a good experimental tool 
for answering the fundamental questions regarding the 
effect of electron-electron collisions. Indeed refined mea- 
surements of thermopower of short ID quantum wires 
have been performed yielding a reasonably good agree- 
ment with the free-electron theory . 27 ' 28 ' 29 It remains an 
open question whether the accuracy of thermopower mea- 
surements is high enough to see the interaction effects in 
longer wires. 
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first order in the transition rate, but here we show it to 
all orders in the interaction. 

We show it explicitly in the Boltzmann equation frame- 
work, however, suspect it to be a general feature of meso- 
scopic systems. Intuitively, the statement means that it 
is the number of particles that passes though the meso- 
scopic system that matters and not their velocity. In con- 
trast to this is e.g. a long ID wire or a bulk metal, where 
a velocity change of the particles is enough to change the 
current. 

To show the above statement explicitly, we formally 
rewrite the Boltzmann equation ([8]) including the bound- 
ary conditions Eq.® as 



f k (x)=f°(e k )+ [ X dx' IkAf] 



for k > 0, (Al) 



la v k 

fk(x) = ft(e k ) + [ dx'^^1 for k < 0. (A2) 

Note that this is not a closed solution of the Boltzmann 
equation, since the distribution function is still contained 
inside the collision integral. However, this rewriting en- 
ables us to find the current without finding the distribu- 
tion function first, i.e. by inserting Eq. (|A1 |A2[) into the 
current definition 



I = 



(~e) 
L 



^2v k f k (x) 



(A3) 



ak 



and obtain (after a few manipulations) : 



a k >0 



=0 



(~e) 
L 



fax ^2 Zk*if\ + { ~r r^'E^'W 



EE 7 (0) + I (int) 



(A4) 



APPENDIX A: SCATTERING PROCESSES 
CONTRIBUTING TO THE CURRENT 

In this Appendix, we show that the particle current 
change due to electronic scattering if and only if the scat- 
tering changes the number of left and right moving elec- 
trons. In the main text (see Eq. (|22|) ). this was shown to 



where the ^-dependent part can be seen to be zero by 
changing variables. We note the cancellation of the ve- 
locity in the distribution function Eqs. l|Al|) . l|A2|l and 
the current definition Eq. 1|A3(1 . which is the origin of 
the statement we are showing (as in the first order calcu- 
lation). A similar cancellation occurs in the Landauer 
formula thus relating the transmission to the conduc- 
tance. By using the explicit form of the collision integral 
Eq. iflOj) the current from the interactions is 



= LA f L dx J2 E W W2 >v [A/2/3(l - - /2')(1 - h>) - fvfvMl - - /a)(l - h)} ■ 

JO (7x0-2(73 kl<0,k 2 ,k 3 

°V°VoV kl ,k 2 ,k 3 , 

(A5) 

We can divide the summation over k quantum number into positive and negative intervals as in the main text (see 
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section HlB[) , The essential point is now, that all terms that have the same number of positive (and negative) intervals 
for the primed and unprimed wave numbers k are zero. In other words, if the number of left and right moving 
electrons does not change then the contribution is zero by symmetry of the transition rate. We show this cancellation 
in practice by an example (using the notation of Eq. (|2"Tj) ): 



E E [/i/ 2 / 3 (l - A<)(1 - /a')(l - h>) - A'/2'/3'(l - A)(l - /a)(l - fa)] (A6) 

spin — I — 

= EE ^23;l'2'3' [/l/ 2 / 3 (l - fv)(l ~ f V ){l ~ h>) ~ fvfvfv<X ~ - / 2 )(1 - /a)] 

spin — + — 

= E E W / 123;l'2'3'/l/ 2 / 3 (l - - /2')(1 - /aO "EE Wka3;l'2'3'/l'/2'/3'(l - /l)(l - / 3 )(l - / 3 ) = 0, 



spin — h— spin — h — 

— | — — | — interchange (123)*=>(1'2'3') 



interchanging 1' and 2' at the first equality using 
Wi2 3; i'2' 3 ' = Wi23;2'i'3' and interchanging (123) <-> 
(1'2'3') in the second term as indicated. Thereby we have 
shown to all orders that to change the current by elec- 
tronic interactions the number of left and right movers 
have to change. 

The statement is not limited to only three particle scat- 
tering and can be show equivalently for pair interaction 
including several bands, electron phonon coupling or any 
other interaction with the same kind of symmetry under 
particle interchange. Furthermore, the statement is still 
true if the collision is non-local in space, since that only 
introduce some spatial integrals in the collision integral 
that can be handled similarly. Note however that the 
distribution function can be changed by processes that 
does not change the number of left and right movers. 



APPENDIX B: DETAILED CALCULATION OF 
THE THERMOPOWER AND CONDUCTANCE 
CORRECTION DUE TO THE THREE PARTICLE 
SCATTERING 

The purpose of this Appendix is to calculate / lnt in 
Eq. {39D 



7 int ~3(-e)^ A i23;i'2' 3 < 

spii 

AT 



spin ++- 



k B T* 



{-£3 + £ 2 ' + £3' - £ F } 



eV 



(Bl) 



in the low-temperature limit, T <C T F , step by step to 
find the prefactor given in Eq. ((5TJ) . As already men- 
tioned, we preform the calculation with the scattering 
process seen in Fig. [2][a) in mind. Therefore we use the 
Fermi functions as given in Eq. I|40p and the matrix el- 
ement entering in the scattering rate from Eq. lj4lj) and 
Eq. l(42|) . i.e. using a quadratic dispersion. 

We preform the summation over all the k in Eq. (|Blj) 
in the following way: First of all, we note that due to the 
momentum and energy conservation in the interaction 



process described, the scattering of fe 3 to fc 3 / has to be 
from above to below the Fermi level, i.e. 



k 3 < 



-k F < k 3 < 



~k F -k 3 )0(k F + k 3/ ). (B2) 



This is due to the signs of ki and and can be under- 
stood as a sign of the difference between the curvature of 
the dispersion near the bottom of the band and near the 
Fermi level. Next, we introduce the momentum transfer 
around the Fermi level qi = k^ — ki for i = 1,3 and using 
the momentum conservation to do the ky summation, we 
obtain 



H 



E E(-) 



(B3) 



fci>0,fc 2 >0,fc 3 <0 qim 



remembering the constraint ky = k\ + q\ > 0, ki' = 
^2 + 51+93 < and k 3 i = k 3 + q 3 < 0. The Fermi 
factors fi (1 — f v ) and /3 1 (1 — f 3 ,) restrict the momentum 
transfer q\ and q 3 to be much smaller then fc F and the ki 
and k 3 to be near the Fermi level for the process in mind. 
Therefore we can use the Fermi functions /f (1 — /{/) to 
do the summation over k\. Assuming slow variation of 
the scattering rate over a range of q\ -C k F at the Fermi 
level, the k\ summation becomes 



J2 0(k F - h)9(h + qi - k F ) = (B4) 

fcl>0 

Similarly the k 3 summation is done using the phase space 
constraint in Eq. (|B2|) 

J2 0(-*p - k 3 )9(k F + ky) = ^q 3 8(q 3 ). (B5) 

fe3<0 

We see that since k\ and k 3 are restricted to the Fermi 
level, we can insert k\ ~ k F and k 3 ~ — fc F in the rest of 
the integrand. To do the k^ > summation, we use the 
energy conservation contained in the scattering rate. It 
is rewritten as (inserting k\ = k F and k 3 = — k F ): 

r , ml 
5{ey + e 2 > + e 3 ' - ex -82- £3) 



x S 



gl - g3 _ 1 

9i + <?3 2 



(91 + 93) 



H 2 \qi + 



13\ 



2 qi + q3 



(B6) 
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We have now done the summation over ki, k 2 and k 3 
and are left with the summation over qi and q 3 of the 
scattering rate, some Fermi functions and the phase fac- 
tors described above. To this end, we introduce w(9i,93) 
by inserting k\ = k F and £3 = — £; F in Eq. l(42|) 



V(ll',22',33') 



^a-^i ,er, $0 



AL 2 h 2 /i 



u( qi ,q 3 ) (B7) 



for a general symmetrized interaction V q = V q + V- q . 

Furthermore, we collect the exponential tales of the 
Fermi functions Eq. l|40bM0cp . 



(1-/^/3(1-/3") = e 
defining 

£(91,93) = £2' ~£3 + £3> 
2 



(e 2 i -S3+E31 -£ F )/k B T 



(B8) 



(B9) 



1 



-hv F ( qi + q 3 ) 



-hvi 



<h ~ 93 

q 3 + qj 2 
(91 -£i)(?l + 9 3 ) , fl 2 ql(2q 2 1 +2q 1 q 3 + ql) 



(«3 + 9l) 5 



2m 



(9i + 93) 2 



inserting ky — k 2 — 9i — 93, from the energy conser- 
vation Eq. (|B6|) and fci = fc F and £3 = — k F . Therefore 
we finally get the interacting contribution to the current 
in Eq. ((HI]) as: 



r int _3(-e)Lm 3 



OO POO 

dgi/ dg 3 
Jo 9i 



9i 93 



K<Zi,9 3 )| 2 (BIO) 



32ir 4 h 7 J J q x + q 3 
X 8(k F - k 2 )6(k 2 )9(-k 2 q 3 )6(k F - q 3 ) 



x e 



(«(9l,93)-E F )/fc B r 



AT 

k^F 2 



(£(91,93) 



k^f 



Here only the step functions, that restricts the integral 
are included. Next we introduce the dimensionless inte- 
gration variables Qi = qi/k F for i = 1,3 and the dimen- 
sionless functions 



S(Qi,Q 3 ) 
in the integral: 
I 



U(Qi, Q 3 ) = k F u(k F Qu k F Q 3 ), 
£(fc F Qi, k F Q 3 ) 



- 3 ^ Lm " e^U<lQ 3: QlQa 



327T 4 fi 7 fc B 



Qi+Q; 



(Bll) 
(B12) 

|c/(Qi,Q 3 )l 2 



xe T 



-^-(S(Q 1 ,Q 3 )-1) 



(B13) 



where .4 is the integration area shown in Fig. Ufa). Note 
that this expression is valid for a general interaction V q 
and that it is not possible to extract a power law in tem- 
perature times some integral by defining new integration 
variables. 

To proceed, we consider the low temperature limit 
T/Tp <C 1 by using the method of steepest decent: Due 




FIG. 4: (a) The integration region A for the integral Eq. (|B13I) 
to calculate the current due to interactions. The two 
boundaries for the integration area close to the origin stem- 
ming from the signs of ki and ky are indicated, (b) The 
H(Qi, Q3) = £(fcpQi, k F Qz)/e F function, a dimensionless ver- 
sion of £(qi, 53) Eq. l|B9p . important in the calculation using 
the method of steepest decent. 



to the exponential function e^""^ 1 '^ 3 ' the maximum 
of 3{Qi,Q 3 ) will dominate the integral for T/Tp — > 0, 
since E(Qi,Q 3 ) < 0. The maxima are 5(0,0) = 
and S(0, 1) = and E(Qi 7 Q 3 ) is shown in Fig. [D[b). 
For a decreasing interaction the area of Q\ <C 1 and 
Q 3 <C 1 dominates even though the integrand is zero for 
Qi = Q3 —* 0. Therefore we expand around the maxi- 
mum (<3i,Q3) = (0,0) to get the lowest order result in 
T/Tp. In view of the integration region Fig. H)Ja), we 
use Q 3 = Qi = Q in the integral Eq. l|B13j) and thereby 
do the Q 3 integral using the approximate limits seen in 
Fig. Ufa), i.e. 



Q1+3Q1 

WQ 3 = 4Ql 



(B14) 



To model the symmetrized potential V q for small q, we 
include the deviation from a constant as described in 
Eq. {48]). This gives 



Qi+Q: 



\U(QuQ 3 



Q 3 =Q!=Q 



1^,9 

9o 



(B15) 
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to lowest order in Q. In the exponential we keep H to 
lowest order in Q, i.e. 



to lowest order in temperature. This is the result stated 
in the text in Eq. (j5lj) . 



e ^H( Qll Q 3 ) ^ e -W. (B16) 

So using the lowest order in Q in the integrand (leading 
to lowest order in T) the interacting contribution to the 
current is: 



int = 3(-e)Lm 3 

327T 4 fi 7 ^ 



x 4Q 2 e- 2< ^ 
8505 



dQV Q 4 
T F AT 



_2ke(^ fc F ) 4 



(o-i) 



M 9 



'2048tt 4 



T 





>'(£ 


7 


AT eV~ 


Wo, 




) 





(B17) 
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